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Abstract 


Nomenclature 


A new procedure seeks to combine the thin- 
layer Navier-Stokes solver LAURA with the parabolized 
Navier-Stokes solver UPS for the aerothermodynamic so- 
lution of chemically-reacting air flowHelds. The interface 
protocol is presented and the method is applied to two 
slender, blunted shapes. Both axisymmetric and three- 
dimensional solutions are included with surface pres- 
sure and heat transfer comparisons between the present 
method and previously published results. The case of 
Mach 25 flow over an axisymmetric six degree sphere- 
cone with a non-catalytic wall is considered to 100 nose 
radii. A stability bound on the marching step size was 
observed with this case and is attributed to chemistry 
effects resulting from the non-catalytic wall boundary 
condition. A second case with Mach 28 flow over a 
sphere-cone-cylinder-flare configuration is computed at 
both two and five degree angles of attack with a fully- 
catalytic wall. Surface pressures are seen to be within 
five percent with the present method compared to the 
baseline LAURA solution and heat transfers are within 
10 percent. The effect of grid resolution is investigated 
in both the radial and streamwise directions. The pro- 
cedure demonstrates significant, order of magnitude re- 
ductions in solution time and required memory for the 
three-dimensional case in comparison to an all thin-layer 
Navier-Stokes solution. 


* Aerospace Technologist, Aerothermodynamics Branch, Gas 
Dynamics Division. Member AIAA. 

t Associate Professor, Department of Aeronautics and Astro- 
nautics. Member AIAA. 

Copyright ©1995 by the American Institute of Aeronautics 
and Astronautics, Inc. No copyright is asserted in the United 
States under Title 17, U.S. Code. The U.S. Government has 
a royalty-free license to exercise all rights under the copyright 
claimed herein for government purposes. All other rights are re- 
served by the copyright owner. 


C 

c p 

D a 

E t 

h 

ho 

H 

k 

kb 

k f 

h 

^ eq 

L e 

M 

M 

n 

N 

P 

q 

R 

Rn 

Re 

s 

T 

U, V, w 

X, Y, Z 
Y b 

V shock 

a 

8 b 

P 

C y, C 

p 


mass fraction 

frozen specific heat at constant pressure 
Damkohler number 
total energy 
enthalpy 

reference enthalpy at the reference temperature 

total enthalpy 

conductivity 

backward reaction rate constant 

forward reaction rate constant 

equilibrium constant 

Lewis number 

Mach number 

molecular weight 

normal distance from the wall 

number of moles 

pressure 

heat transfer rate 

universal gas constant, 8314.3 J/kg-mole-K 
nose radius 
Reynolds number 

distance measured along body surface 
temperature 

cartesian velocity components 

cartesian coordinates 

body location 

shock location 

angle of attack 

body angle 

viscosity 

curvilinear coordinates 
density 

characteristic time 


Subscripts: 
i LAURA 


species 
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wa ii value at the wall 

oo freestream value 


Introduction 

International interest in a space station, the possi- 
bilities for human exploration to other planets, and the 
advancing age of the space shuttle fleet have all brought 
the issue of advanced launch and reentry vehicles to the 
forefront. A critical design point for these vehicles is 
during hypersonic reentry, when peak heating rates oc- 
cur and aerodynamic control effectiveness may be altered 
due to flowHeld phenomenon unique to the high-altitude, 
high-velocity conditions. The high temperatures and 
high convective velocities relative to reaction times cre- 
ate an environment where chemical nonequilibrium ef- 
fects can be significant. Accurate aerothermodynamic 
predictions during this part of the reentry trajectory are 
essential for sizing both the thermal protection system 
and aerodynamic control surfaces. Ground based tests 
simulating these flight conditions, including considerable 
nonequilibrium effects, are difficult to perform. Flight 
tests can be prohibitively expensive. 

Two popular computational approaches for obtain- 
ing aerothermodynamic predictions on these classes of 
vehicles are to solve the thin-layer Navier-Stokes (TLNS) 
equations or the parabolized Navier-Stokes (PNS) equa- 
tions. TLNS is derived from the full Navier-Stokes equa- 
tions by neglecting viscous terms in the streamwise and 
crossflow directions. The assumptions inherent in the 
TLNS equations are often acceptable for a wide class 
of conditions and configurations, including cases of hy- 
personic, chemically reacting flow. Excessive compu- 
tational requirements can become a drawback to using 
TLNS as the entire solution domain is relaxed in time. 
Complex configurations 1 can tax computer memory with 
millions of grid points, and solution times may be mea- 
sured in CPU days. In addition, solving for chemical 
nonequilibrium can, for some algorithms performing ex- 
act matrix inversions, increase the computer memory 
and time requirements by the cube of the number of 
species considered . 2 

The PNS equations are obtained from the full 
Navier-Stokes equations by neglecting the time deriva- 
tives and the streamwise viscous derivatives. Limited to 
flowHelds with streamwise supersonic flow outside the 
boundary layer, no streamwise separation, and weak 
streamwise pressure gradients in the subsonic region, 
PNS algorithms are well suited for solving sharp-nosed, 
slender-body supersonic/hypersonic configurations. Be- 
ing space marching and steady state, PNS formulations 
can realize appreciable decreases in both computational 


time and memory requirements relative to TLNS algo- 
rithms. The principle difficulty in applying the PNS 
equations to the class of vehicles considered here is 
that commonly the algorithms cannot solve blunt-body 
flowHelds, and most reentry vehicle designs incorporate 
blunted nose and leading edge regions in order to reduce 
peak heating rates. 

The present study looks to combine two well- 
established computational codes, one a TLNS algorithm 
and the other PNS, for the solution of chemically- 
reacting, hypersonic flowHelds. The technique is suc- 
cessfully applied to blunted, multi-conic geometries at 
both zero and non-zero angles of attack. Different sets 
of freestream conditions are considered, and the effects 
of wall catalycity are investigated. Challenges and ob- 
stacles to the consistent integration of the two codes are 
observed and comments regarding the applicability and 
limitations of the procedure are documented. 

Prelude to the Present Method 

Recently, efforts have been made to combine the 
TLNS and PNS approaches in order to get timely, ac- 
curate hypersonic viscous solutions while circumventing 
some of the above mentioned limitations. Weilmuenster 
and Gnoffo 3 proposed a multi-block solution procedure, 
in which the domain is divided into blocks ordered in the 
streamwise direction. The general idea is to march these 
blocks downstream, analogous to the PNS approach of 
marching two-dimensional planes, and to solve the inte- 
rior of each block with TLNS. This procedure principally 
attacks the memory requirements inherent in obtaining 
a full-body TLNS solution by splitting the domain, but 
does not decrease the time required to obtain the solu- 
tion since TLNS remains the governing equations. 

The TLNS code used by Weilmuenster is the 
Langley Aerothermodynamic Upwind Relaxation Algo- 
rithm (LAURA ). 4-7 LAURA is a finite- volume, shock- 
capturing, hyperbolic equation solver with second-order 
spatial accuracy for the steady-state solution of viscous 
or inviscid hypersonic flows. The scheme employs a 
point implicit relaxation strategy with the upwind flux- 
difference splitting of Roe . 8 The right-hand-side (RHS) 
of the equations are formulated according to Yee 9 with 
the entropy condition of Harten . 10 Perfect gas, equilib- 
rium air, and nonequilibrium air calculations can all be 
performed. 

Greene 11 has extended the LAURA code into a PNS 
version. With this method, LAURA-TLNS is used on 
the blunt-nose portion of a hypersonic vehicle. At a point 
in the flowfield consistent with the PNS equations, the 
transfer is made to LAURA-PNS, which is then marched 
down the remainder of a slender vehicle afterbody. This 
particular formulation, being a TLNS extension, is lo- 
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cally iterative in pseudo-time steps, and its performance 
suffers from arriving as a PNS solver via a TLNS algo- 
rithm, rather than being a code that was optimized as 
a PNS solver from inception. Thus, while this method 
significantly reduced the memory required to obtain a 
solution, it was not able to reduce solution time to the 
level desired. 

Upwind Parabolized Navier-Stokes Solver 
(UPS) 12-20 is an upwind, finite- volume, state-of-the-art 
PNS code with chemical nonequilibrium capability. It 
is second-order accurate in the crossflow plane and first 
order accurate in the marching direction. The equations 
are approximately factored and solved implicitly, with 
the approach of Vigneron et al , 21 employed to suppress 
departure solutions. 

UPS was identified as a code that, when com- 
bined with LAURA, might provide the tremendous re- 
duction in vehicle solution time originally sought with 
the LAURA-TLNS/LAURA-PNS method. The present 
method seeks to combine LAURA and UPS for a consis- 
tent solution procedure for air flows in chemical nonequi- 
librium. 

Previously, UPS has been joined with the TLNS 
code CNS by Lawrence et al. 15 for perfect gas computa- 
tions. Nonequilibrium solutions are presented by Buelow 
et al. 22 and Muramoto 23 using UPS with the TLNS code 
TUFF. LAURA has the advantage over TUFF in that it 
can handle generic, three-dimensional geometric shapes, 
as are encountered with real vehicle configurations, and 
is an upwind, finite volume method, like UPS. 


Present Method 

A combined LAURA-UPS solution procedure has 
been implemented by Wood and Thompson 24 for per- 
fect gas and equilibrium air flows. That study included 
detailed solutions for an axisymmetric perfect gas case 
and a three-dimensional equilibrium air solution for the 
Reentry F vehicle, 25 including turbulence. Generally 
good results were seen with the combined method, and a 
very significant reduction in solution time was achieved. 
The extension of this procedure to nonequilibrium air 
calculations, however, is not straightforward, because 
while both UPS and LAURA use the same equilibrium 
air curve fits, they do not use the same chemistry models 
for nonequilibrium air. 

The TLNS LAURA solutions were sought using 
a chemical nonequilibrium, thermal equilibrium kinetic 
model with a seven-species air model. The PNS solu- 
tions were obtained with UPS using the seven species, 
seven reaction nonequilibrium air chemistry model with 
the reaction rates of Blottner et al. 28 


Modifications for Compatibility 

Changes made to the LAURA pre-processor for 
compatibility with UPS focus mostly on grid genera- 
tion. The grid on the cone portion of a sphere-cone was 
changed from being body normal to being axis normal 
so as to facilitate space marching on slender bodies. The 
spacing normal to the body in the initial grid was mod- 
ified so as to better capture the bow shock for vehicles 
with very slender afterbodies. The number of cells solved 
on spherical nosecaps was reduced to 12. 

The wall boundary conditions in LAURA were 
changed to correspond with the UPS wall boundary con- 
ditions by switching from the standard LAURA bound- 
ary conditions to the primary alternate boundary condi- 
tions. The standard LAURA viscous wall boundary con- 
ditions apply the wall values, i.e., zero velocity, fixed wall 
temperature, etc., at the center of an image cell below 
the vehicle surface. The UPS approach is to use reflected 
boundary conditions for the image cell, so as to apply the 
boundary conditions to be at the cell face defining the 
wall. The UPS approach is considered to be a higher- 
order method than the default LAURA boundary impo- 
sition. However, the LAURA default boundary condi- 
tions were found to be more robust than the reflected 
boundary conditions, so the LAURA solutions were first 
partially converged with the standard boundary condi- 
tions, and then switched to the reflective boundary con- 
ditions during the later stages of convergence after the 
flowHeld had stabilized. This switch is usually made at 
the same time spatial second-order accuracy is enforced. 

Of the five kinetic models available in LAURA, the 
15 reaction model of Kang et al. 29 was chosen as being 
the closest match with the rates of Blottner 28 in UPS. 
Two further parameters were toggled from the default in 
LAURA to better deal with slender-body configurations. 
The eigenvalue limiter was set to be scaled by the cell 
aspect ratio and the upwinding of the surface properties 
was turned off. 

The principle change made to UPS involved the 
restart Hie. A jump in properties was observed during 
nonequilibrium restarts. This was tracked to the use of 
a freestream value of the mixture molecular weight when 
initially decoding the temperature from the energy and 
species concentrations, prior to marching. The remedy 
was to read the local mixture molecular weight into the 
standard restart file. 

Remaining Differences Between the Codes 

Some differences in the chemistry models remain 
between LAURA and UPS. Algorithmically, LAURA 
solves the chemistry equations with a fully-coupled pro- 
cedure while UPS uses a loosely-coupled approach, but 
with the option for local subiterations to get a close ap- 
proximation to a fully-coupled scheme. The two codes 
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Fig. 1 Comparison of molecular oxygen species 
viscosities. 

compute the species enthalpies with fundamentally dif- 
ferent approaches, as LAURA uses polynomial curve fits 
while UPS uses interpolated table look-ups. This pre- 
vents the exact matching of species concentrations, inter- 
nal energy, and temperature between the codes, though 
the magnitude of the difference is considered to be small 
enough to not prohibit the interfacing of the codes. 

Further differences exist in the way each code com- 
putes the bulk thermodynamic and transport properties. 
This leads to small mismatches between the codes for pa- 
rameters such as viscosity and speed of sound. One ques- 
tion this raises is whether to match the non-dimensional 
freestream quantities Mach number and Reynolds num- 
ber between the codes, or to match the dimensional 
freestream velocity and density. For high Reynolds num- 
ber, hypersonic applications where Mach number inde- 
pendence has been reached, the decision made here is 
to match the dimensional freestream conditions. As ex- 
amples of the differences in the transport property com- 
putations, plots of viscosity versus temperature are pre- 
sented for molecular oxygen, Fig. 1, and nitric oxide, 
Fig. 2. The computations of both LAURA and UPS are 
presented along with the recommendations of Gupta et 
al., 30 who conducted one of the most recent studies into 
transport property computations. Generally, the other 
species viscosities match fairly well over the tempera- 
ture range 1000-30,000 Iv, with the molecular nitrogen 
viscosity computations presented in Fig. 3 as a typical 
example. Sample computations of thermal conductivi- 
ties performed for typical near wall conditions resulted 
in a 4-5 percent higher value from UPS than LAURA. 
R is difficult to predict a priori what effect these dif- 
ferences would have relative to solutions from the two 
codes. 



Fig. 2 Comparison of nitric oxide species viscosities. 



Fig. 3 Comparison of molecular nitrogen species 
viscosities. 


Interface Protocol 

The interface procedure between LAURA and UPS 
begins with the standard LAURA restart file for a 
converged chemical-nonequilibrium solution. From the 
LAURA restart file, a crossflow data plane is extracted to 
become the UPS starting plane. Currently, this plane is 
chosen at least three cells upstream of the final LAURA 
solution plane in order to avoid possible contamina- 
tion from the extrapolated outflow boundary conditions. 
The variables available in the LAURA restart file are: 
the three velocity components, temperature, the seven 
species densities, and the finite volume grid, 

[u,v,w,T, p s ,x,y,z\ 

The variables needed by UPS to start are: mixture 

density, the three momentum components, total energy, 
mixture molecular weight, species mass fractions, and 


4 



the starting plane in finite volume form, 

\p, pu, pv, pw, E t , M, c s ,x, y, z\ u 

In the equations which follow, a subscripted “1” is used 
to indicate a LAURA variable or quantity, while the sub- 
script “u” refers to the corresponding UPS parameter. 

The variables required by UPS are obtained from 
the LAURA variables in the following manner. The grid 
is transformed according to the transformation of the 
physical coordinates as, 

Xu = -Zi , Vu = XI , z u =yi (1) 

The total density is found from summing the species 
densities, 

Pu = 'y ] Ps,i (2) 


the Vigneron condition’s limitation on the pressure gra- 
dient near the wall set up oscillations that restricted the 
stability of the marching UPS solution. The fix to this 
problem was to pass the LAURA temperatures directly 
through to the UPS species enthalpy interpolated table 
look-ups, then to complete the computation of the total 
energy as described above, 

E t , u = + vf + wf) + H u (7) 

^ Pu 


where, 



and, 


u — Cp.u Pi T h 


<■ 0 ,s ,u 


( 8 ) 

(9) 


The three components of momentum are obtained from 
the velocity components and the total density, 

PUu = Pu • (~wi) , pv u = p u -ui , pw u = Pu-vi (3) 


Species mass fractions are found by dividing the species 
densities by the total density, 


Cs,u 


Ps,l 

pu 


(4) 


The mixture molecular weight is found by applying the 
perfect gas equation of state to the mixture temperature 
and pressure, 


M u 


PuRTi 

Pu 


(5) 


where the mixture pressure was determined from sum- 
ming the species partial pressures, 


Pu 


p^Yl 


s 


c s ,uRTi 

M s 


( 6 ) 


a step consistent with the assumption, common to both 
codes, that the working fluid is a mixture of ideal gases. 

The UPS total energy now remains to be computed. 
Initially, the effort was made to take the temperature and 
species densities from the LAURA solution, pass them 
through the LAURA enthalpy curve fits, add in the ve- 
locity and species property information, and obtain a 
total energy that would be passed directly to UPS. A 
problem was encountered when UPS took this energy 
and decoded temperature and pressure. The differences 
between the UPS and LAURA enthalpy computations 
lead to differences between the decoded UPS tempera- 
tures and pressures and the original LAURA tempera- 
tures and pressures. These variations, in combination 
with a fixed wall-temperature boundary condition and 


Since both LAURA and UPS are finite volume for- 
mulations, the UPS starting-plane grid is taken at a 
streamwise location corresponding to the location of the 
cell-centered LAURA data. Converting the nondimen- 
sionalizations so that the UPS velocities are normal- 
ized by the freestream speed of sound, rather than the 
freestream velocity as is done in LAURA, and performing 
the curvilinear transformation between the two codes, 

Cu = -m , Vu = 0 , Cu =6 (io) 

completes the interface process. 


Results 

The present method is successfully applied to two 
primary configurations and flow conditions. Case 1 is an 
axisymmetric sphere-cone, chosen to correspond with the 
results of Gupta et al. 31 The nose radius is 0.0381 m and 
the body angle is six degrees. The freestream conditions 
are for Mach 25 at an altitude of 53.34 km (175 kft.). The 
wall temperature is held fixed at 1260 K, with a non- 
catalytic chemistry boundary condition. Case 2 is for 
Mach 28 flow over the sphere-cone-cylinder-flare configu- 
ration studied by Bhutta et al. 32 at both two and five de- 
gree angles of attack. This configuration has a 0.1524 m 
spherical nosecap followed by a nine degree cone. Af- 
ter 10 nose radii the cone is followed by a cylinder and 
then a five degree flare, each of 10 nose radii length. 
The freestream conditions correspond to an altitude of 
83.8 km, (275 kft.), at a Reynolds number per meter of 
6148. The wall temperature for this case is 833 K and 
a fully-catalytic boundary condition is employed. Ta- 
ble 1 presents a summary of the nominal conditions for 
the two cases. For all calculations the freestream species 
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Table 1 Nominal conditions for Cases 1 and 2. 



Case 1 

Case 2 

Configuration 

sphere-cone 

blunted multi-conic 


25 

28 

Re (nr -1 ) 

3.95 xlO 5 

6148 

Altitude (km) 

53.34 

83.8 

Rn (m) 

0.0381 

0.1524 

0b (deg) 

6 

9-0-5 (10 R n each) 

T W all (K) 

1260 

833 

Wall catalycity 

none 

fully 

a (deg) 

0 

2, 5 



concentrations were set at, 


cjv 2 


0.767 

c 'o 2 


0.233 

cn 


6.217 x 10 _2 ° 

co 


7.758 x 10 -9 

cno 


4.981 x 10- 5 

c NO + 


4.567 x 10“ 24 


The seventh specie, electrons, are found from a charge 
balance with the ionized nitric oxide, 


M e — Nno+ ( 12 ) 


Case 1 

A viscous, second-order accurate TLNS LAURA so- 
lution was obtained for Case 1 with chemical nonequi- 
librium, thermal equilibrium, and a non-catalytic wall 
condition, implemented in both codes as, 


dc s 

dn 


= 0 

wall 


(13) 


i.e. t the mass fractions of the image cells are set equal 
to the mass fractions of the first cell outside the wall. 
The axisynnnetric LAURA computational grid contains 
64 cells normal to the body and 28 cells in the stream- 
wise direction, extending five nose radii to 0.19 nr. This 
grid was adapted using the standard LAURA grid adap- 
tion routine. Figure 4 displays the final LAURA grid, 
for clarity showing only every fourth point in the body- 
normal direction. For consistency, Fig. 4 and all subse- 
quent figures use the UPS coordinate system. The loca- 
tion win -re the UPS starting plane was extracted from 
the LAURA solution is indicated in Fig. 4. That por- 
tion of the LAURA grid downstream of the UPS starting 
plane was supplied as an external grid to UPS. Since the 
UPS marching step size was smaller than the LAURA 
cell sizes shown in Fig. 4, the LAURA grid was linearly 


Fig. 4 Case 1 LAURA computational grid, showing 
every fourth body-normal point. 


interpolated in the streamwise direction to obtain the 
actual UPS grid. This is the standard UPS approach 
for handling external grids. A significant overlap of the 
solution domains was deliberately chosen for this case 
to allow for a direct code-to-code comparison between 
LAURA and UPS. In general, an overlap of this size is 
not required by the combined procedure. 

The UPS solution was carried out 100 nose radii to 
3.81 nr by extending the external grid downstream in a 
conical extrapolation. The grid was moderately adapted 
to the solution in the body-normal direction as the solu- 
tion proceeded, in such a way as to maintain the origi- 
nal grid spacing at the wall while linearly stretching the 
outer 60 percent of the grid. This adaption routine is 
currently not fully integrated into the version of UPS 
used here, and relies upon the user to provide the neces- 
sary stretching parameters. 

Figure 5 displays Mach contours from the LAURA 
and UPS solutions, covering the overlap region to five 
nose radii. The location of the UPS starting plane is 
indicated, and the UPS Mach contours are overlaid upon 
the LAURA Mach contours downstream of that point. 
Excellent agreement is seen between the present method 
and the LAURA-only solution. 

Figures 6 and 7 plot surface pressures, normalized 
by twice the freestream dynamic pressure, versus the 
streamwise distance measured along the surface, normal- 
ized by the nose radius. The viscous shock-layer (VSL) 
solutions of Gupta 31 are included for comparison. The 
VSL equations employ a further approximation to the 
governing equations beyond the PNS equations to allow 
solution marching in both the streamwise and circumfer- 
ential directions. Figure 6 is a close-up on the interface 
region, extending to 10 nose radii. The Gupt.a-VSL so- 
lution extends the full 10 nose radii, while the LAURA 
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Fig. 5 Mach contours:; tIPS solution overlaid upon 
LAURA solution. 



Fig. 6 Case 1 surface pressures — the interface region. 
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Fig. 7 Case 1 surface pressures to 100 nose radii. 


Fig. 8 Case 1 surface heating — the interface region. 

solution was terminated at six nose radii in this plot. 
The UPS solution was initiated at two nose radii and 
extends to 10 nose radii. Excellent agreement is seen 
between the UPS and LAURA solutions. The Gupt.a- 
VSL solution is seen to agree very well with the LAURA 
and UPS solutions outside of the region of sphere-cone 
tangency, where Gupta-VSL predicts higher pressures, 
probably due to the surface curvature smoothing em- 
ployed in this VSL code. A slight pressure bump ap- 
pears in the UPS solution at five nose radii. The cause 
for this is not known at this time, and may be a residual 
of the LAURA-pPS interface. However, the effect is lo- 
calized and does not appear to influence the downstream 
solution. 

Figure 7 extends the surface pressure plot out to 100 
nose radii, capturing the overexpansion and recompres- 
sion regions. There is a maximum difference between 
the UPS and Gupta-VSL solutions of 3-4 percent in the 
reconrpression region. As in Fig. 6, the LAURA solution 
was terminated at six nose radii. 

Surface heat transfer results for LAURA, UPS, and 
Gupta-VSL are presented in Figs. 8 and 9. Figure 8 
plots the interface: region out to a distance of six nose 
radii. Similar trends are seen in the heating as were seen 
for the pressure in this region. The heating at the in- 
terface between the. LAURA and UPS codes picks up 
smoothly, but there is a bump in the UPS heating be- 
tween four and five nose radii, corresponding to the pres- 
sure bump discussed earlier. The Gupta-VSL heating is 
elevated above the LAURA-UPS heating in the region 
of the sphere-cone juncture. 

Figure 9 carries the present method and Gupta-VSL 
heating out to 100 nose radii. Note that the LAURA 
heating terminates at six nose radii. A noticeable differ- 
ence exists between the UPS and Gupta-VSL solutions 
that persists from the overexpansion region on down- 
stream. The Gupta-VSL results are consistently 18-22 
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Fig. 9 Case 1 surface heating to 100 nose radii. 
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Fig. 10 Case 1 atomic oxygen mass fraction profiles 
at X/R n = 5. 


percent lower than the UPS heating. More investigation 
is required to understand why there is this level of differ- 
ence between the solutions, but, while Figs. 1 and 2 show 
good agreement between the UPS and Gupta viscosities, 
there are differences in other aspects of the kinetic mod- 
els which may be contributing to the heating disparity. 

Looking specifically at reacting chemistry effects, 
Fig. 10 profiles the atomic oxygen mass fraction ver- 
sus normal distance from the surface, as a fraction of 
the shock layer, at an axial position five nose radii from 
the noset.ip. The profiles from the LAURA and UPS so- 
lutions are seen to be similar, with a difference in mass 
fraction at the surface of two percent. The mass fraction 
gradients at the surface are seen to be zero, as defined 
by the non-catalytic wall assumption. 

This difference in oxygen mass fraction at the sur- 
face becomes critical in realizing the difficulties encoun- 
tered in obtaining the combined LAURA-UPS solution 
for this particular configuration. While the Blottner and 


Table 2 Reaction rates for ionized nitric oxide. 


kf 

Kang 1.4 x lO 6 ^- 5 exp-^ 6.7 x W 21 T~ 15 

Blottner 9.03 x 10 9 T°- 5 exp- 22 ^ 1.8 x lO 19 ^ 1 0 


Kang reaction sets are similar or identical for most reac- 
tions, a significant difference in the equilibrium constant 
can occur in the equation controlling production of ion- 
ized nitric oxide, 

N + 0 — NO + +e~ (14) 

Table 2 lists the forward and backward rates for Eqn. 14 
from the two kinetic models. At a temperature of 1280 
Iv, an average temperature for a Case 1 surface cell, the 
Blottner equilibrium constant, 

ke q = ^ (15) 

for this reaction is 2.33 x 10 _lb , while the Kang equi- 
librium constant is two orders of magnitude lower at 
6.58 x 10 -18 . Under the flow conditions for this case, both 
atomic oxygen and atomic nitrogen concentrations at 
the surface are large, with the flow consisting of roughly 
equal parts atomic oxygen, atomic nitrogen, and molecu- 
lar nitrogen near the wall. The net result is that the UPS 
solution produces significantly more ionized nitric oxide 
relative to the starting solution provided by LAURA, 
and at a fast rate. This creates a marching step-size sta- 
bility restriction characterized by the Damkohler num- 
ber, 

D a = Tflow (16) 

T reactions 

see Ref. 33 pp. 149-154, which was found to be exacer- 
bated by a tight grid spacing near the wall. 

A compromise was sought whereby the LAURA grid 
was modified to double the cell size of the first grid cell, 
which sets a nominal cell Reynolds number of two at 
the wall. This was found to still allow accurate reso- 
lution of gradients at the wall while somewhat relaxing 
the Damkohler imposed marching stability restriction. 
In this case the non-linearity inhe®nt. in the chemical 
reactions allows for marching steps more than twice as 
large as were possible with a wall cell Reynolds number 
of one. Larger grid spacings at i lie wall were found to 
be too coarse to provide suitable LAURA solutions. 

The LAURA solution for this case was converged 
through an L'> norm of the residual of seven orders of 
magnitude in 2200 iterations. The total CPU time on a 
Cray 2 was 1411 seconds. Figure 11 contains the con- 
vergence history of the LAURA solution. One caveat to 
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Fig. 11 Case 1 LAURA convergence history. 

this performance is that the solution was begun from a 
converged solution for a similar, but not identical, case. 
The sharp spikes occurring early in Fig. 11 are the re- 
sult of grid adaptations, while the later spikes are due to 
shock ringing. 

The UPS solution was obtained on a Cray YMP 
with a final marching step size of 0.25 mm. This is 
a small step size in relation to other cases which have 
been run with the present method, but is a result of 
the previously mentioned marching stability restrictions. 
Muramoto 23 reports using the same marching step size 
for a Mach 20, seven-degree sphere-cone nonequilibrium 
case, with a modified version of UPS, and Tannehill 
et al , 18 report using a step size of 0.2 mm on an axisyrn- 
rnetric cone. The full UPS solution to 100 nose radii 
required 4198 seconds. 

Attempts at Larger Marching Steps 

While the axisymmetric geometry of Case 1 was able to 
be solved by the present method in a reasonable amount 
of computer time with the small marching step size, 
there is concern that a full-sized, three-dimensional ve- 
hicle might require excessive computational resources if 
conditions wefe such that the stability restriction ob- 
served here applied. Several attempts were made to 
enhance the stability of the UPS marching solution 
for the Case 1 conditions. Local chemistry iterations 
were added, second- and fourth-order subsonic smooth- 
ing terms were turned on, the safety factor applied to 
the Vigneron condition was adjusted, and the eigenvalue 
stability parameters EPS A and EPSS in the UPS in- 
put file were changed. Some small stability improvement 
was found by increasing the values of the second-order 
implicit smoothing term, the Vigneron safety factor, and 
the stability parameter EPS A, but not enough to allow 
order-of-magnitude larger step sizes. Direct reprogram- 
ming of either code’s complete chemistry package to ex- 


actly match the other was not attempted, being beyond 
the scope of the present study. 

The UPS solution instability for this case is typically 
manifested by a divergence of the Sell temperatures at 
the wall. It was thought that the reflected boundary con- 
ditions, where the wall temperature is enforced only as 
the geometric average of the image cell temperature and 
the temperature of the first cell above the wall, might 
be contributing to the instability because the wall tem- 
perature is not explicitly enforced. The UPS boundary 
conditions were altered to apply the wall boundary con- 
ditions of no slip, no penetration, and fixed wall temper- 
ature at the image cell center, and a new solution was 
obtained with both LAURA and UPS using this bound- 
ary condition, but no appreciable improvement in the 
stability of the present method was observed. An ef- 
fort to enforce a limiter on the Newton iteration used 
to decode the temperature and pressure given the to- 
tal energy, mixture density, and species concentrations 
also failed to produce a useful relaxation of the stability 
restriction on the marching step size. 

Some attempts at solution smoothing and solution 
modulation were tried with the present method. Sev- 
eral approaches were attempted, beginning by trying to 
march the UPS solution one step, modifying the origi- 
nal interface plane with an under relaxation scheme by 
adding some fraction of the difference between the ini- 
tial starting plane data and the first step solution, and 
repeating in a locally iterative procedure. The idea was 
to allow the solution to relax without creating excessive 
transients. The next attempt tried to march the UPS 
solution while modulating it with the LAURA solution, 
so that the first step was 10 percent UPS and 90 per- 
cent. LAURA, the second step 20 percent UPS and 80 
percent LAURA, and so on. While these attempts had 
some small success in delaying or postponing the insta- 
bility with large step sizes, they were unable to suppress 
the instability enough to solve a significant portion of 
the geometry with large marching steps. More exotic so- 
lution modulation methods were tried whereby the UPS 
domain was split to allow the inviscid, viscous, and near- 
wall regions to relax from the LAURA solution at differ- 
ent rates, but the result was still the same — the marching 
step-size was limited to the millimeter range or less. 

An attempt at a solution was made using the Park 33 
kinetic model in LAURA, with no more justification than 
that it is a readily available option. Perhaps predictably, 
this did not produce any improvement in stability. The 
location of the interface point was varied as well, without 
producing a change in the behavior of the solution with 
the present method. 

Changes to the grid included trying 40, 64, and 128 
points in the body normal direction with nominal cell 
Reynolds numbers at the wall of 0.5, 1, 2, 5, and 10. 
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Fig. 12 Case 1 heat transfers for a fully-catalytic wall. Fig. 13 LAURA Case 2 symmetry plane grid, showing 

every eighth radial point. 


The number of points did not seem to alter the solution 
appreciably for this configuration, but as discussed ear- 
lier the cell size at the wall proved to be very important. 
The tradeoff had to be made between a tight clustering 
at the wall for good gradient resolution and a more rea- 
sonable cell aspect ratio to allow feasible marching step 
sizes. 

A final parametric on the basic Case 1 solution was 
performed by employing a fully-catalytic wall instead of 
the non-catalytic boundary condition. The surface heat- 
transfer results for this case are presented in Fig. 12. 
For this solution a march larger step size was possible 
with UPS, because the fully-catalytic wall condition cre- 
ates a different gas composition in the near-wall region 
which does not involve the ionized nitric oxide reaction, 
Eqn. 14, to the same degree as the non-catalytic solu- 
tion. However, as can be seen in Fig. 12 the heating 
from the present method immediately downstream of the 
interface region does not look good, although the heat 
transfers agree well from 15-100 nose radii with the re- 
sults of Gupta for the same configuration. Interestingly, 
the same behavior in the UPS heating near the interface 
point is reported by Muramoto in Fig. 11 of Ref. 23 for 
an axisymmetric, seven degree sphere-cone with a fully- 
catalytic boundary condition. In discussing his result, 
Muramoto further cites Buelow 22 as another investiga- 
tor who has seen a similar heating behavior with UPS. 


Case 2 

Three-dimensional Hnite-ratfi. chemistry solutions 
were sought for the Case 2 configuration for Mach 28 at 
two and five degree angles of attack. The wall boundary 
condition was set to be fully catalytic, which in both the 
LAURA and UPS versions employed here sets the wall 


species concentrations equal to their freestream values, 


Cs,wall — Cs,oo 


(17) 


Since the species concentration gradients are no longer 
zero at the wall for this case, as they were for the non- 
catalytic solutions, a computation of the diffusive heat- 
ing rate was added to the UPS surface property output 
routine as, 


_ k L e \ u dc s 
~ rt. 2^ s 7h7 


Computations of the diffusive heating for the cases con- 
sidered in the present study showed its contribution to 
the total heat transfer to be a very small percentage. In 
the calculations of Ref. 32 a variable wall temperature 
was employed, but for the present calculations a fixed 
wall temperature of 833 Iv was used. This was chosen as 
a rough average to use for comparison with the results 
of Bhutta. 


Two Degrees Angle of Attack 

The LAURA symmetry plane grid for this case is dis- 
played in Fig. 13, for clarity showing only every eighth 
point in the radial direction. The full LAURA grid con- 
tains 51 streamwise cells, 18 circumferential cells, and 
128 radial cells. The UPS starting plane was extracted 
from the fifteenth streamwise cell in the LAURA solu- 
tion. 

Figure 14 contains both windside and leeside surface 
pressures, normalized by twice the freestream dynamic 
pressure, versus axial distance, normalized by the nose 
radius, for both the full-body LAURA solution and the 
coupled LAURA-UPS solution of the present method. 
The agreement is very good, with the most noticeable 
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Fig. 14 Case 2 surface pressures, windside and leeside. 

difference occurring at the cylinder-flare junction. The 
pressure jump at the flare is much more sharply defined 
with the UPS solution, whereas LAURA predicts a less 
abrupt pressure change. Part of this difference is at- 
tributed to the prevention of upstream propagation of 
pressure waves in the subsonic portion of the boundary 
layer by the PNS code. Also contributing is a somewhat 
coarse LAURA grid in the streamwise direction at this 
point. The UPS solution was marched at a step size of 
0.01 nr, which is about, one-tenth the streamwise length 
of the corresponding LAURA cells at the cylinder-flare 
junction. Remember though that LAURA is second- 
order accurate in the streamwise direction, while UPS 
is only a first-order algorithm in the marching direction. 
R. can be seen that with the present method UPS picks 
up the pressure accurately from the LAURA solution at 
the interface region, located at two nose radii. 

A mirrored pressure contour plot is presented in 
Fig. 15. The left half of the solution is from LAURA 
while the right side is the UPS solution. Both solutions 
are taken from a cross section at 29 nose radii. It can be 
seen that the UPS bow shock is crisper than the LAURA 
bow shock. This feature holds true in general, and is a 
result of the different types and levels of numerical dis- 
sipation used in each scheme. Reasonable agreement is 
seen between the two solutions, 27 nose radii downstream 
of the interface point. 

Axial surface heat transfers are plotted in Figs. 16 
and 17. Along with the results from the two codes in the 
present method, heat transfer results from Ref. 32 for a 
nonequilibrium PNS code, PNSNQ3D, and a nonequilib- 
riurn VSL solver, VSLNEQ, are presented for compari- 
son. In Fig. 16 very good agreement is seen with the 
present method, as the distribution in windside heat- 
ing spans 5-10 percent between LAURA, UPS, and 
PNSNQ3D over the vehicle body. The VSLNEQ results 



Fig. 15 Case 2 pressure contours in the crossflow 
plane at X/R n = 29. 



Fig. 16 Case 2 windside heat-transfer rates. 

are as much as 35 percent lower than the other solutions 
on the cylinder. Looking at the interface region, the UPS 
heating is seen to pick up very well from the LAURA 
starting solution. At the juncture between the cylinder 
and the flare, the UPS solution is seen to capture a more 
abrupt change in heating than the LAURA solution. As 
was the case with the surface pressure, the. cause of this 
difference is attributed to the suppression of upstream 
information propagation by the space marching scheme 
and axial smearing by the LAURA grid. 

The corresponding leeside heat transfers are plot- 
ted in Fig. 17. Leeside heating for this case, with a two 
degree angle of attack, is 40 percent lower than the wind- 
side heating. The same trends between the four solutions 
are seen on the leeside as on the windside, with even 
slightly better agreement. The LAURA and UPS solu- 
tions on the leeside agree to within nine percent. Again, 
the VSLNEQ results are lower than the other heating 
rates. 
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Fig. 17 Case 2 leeside heat-transfer rates. 



Fig. 18 Case 2 LAURA convergence history. 
Timing 

It was seen earlier that for the non-catalytic wall con- 
ditions in Case 1 the present method was limited in 
its computational advantage over the full TLNS solu- 
tion by a marching stability step-size restriction. This 
is not the case for the fully-catalytic surface of Case 2. 
Since the species concentrations are forced to return to 
freestream values at the wall, there is very little atomic 
oxygen and atomic nitrogen at the wall, and hence the 
reaction controlling production of ionized nitric oxide, 
Eqn. 14, is not the factor it was in Case 1. Much larger 
marching step sizes were able to be taken for the sphere- 
cone-cylinder-flare configuration than for Case 1, and a 
substantial reduction in solution time was achieved with 
the present method over a full TLNS solution. Both the 
LAURA and UPS solutions for this case were obtained 
on a Cray YMP. Figure 18 tracks the convergence his- 
tory of the LAURA solution for this case. In this plot 
the residual starts out small and then jumps up abruptly. 
This is part of the initialization and restart procedure, 


and does not represent a converged solution. As with 
Fig. 11, the spikes in the convergence history during the 
first hour are the result of grid adaptations. The later 
spikes are associated with the multi-tasked restart proce- 
dure in LAURA. The total LAURA solution CPU time 
was 4.73 hours, requiring 25 megawords of memory. 

With a marching step size of 0.01 nr, two orders of 
magnitude larger than were possible for Case 1, the UPS 
portion of the solution was obtained in only 776 CPU 
seconds, and required only 2.15 megawords of mem- 
ory. This represents an order of magnitude reduction in 
both time and memory over the full-body LAURA so- 
lution. Results presented in the next section show that, 
with a slight reduction in solution resolution, the present 
method can achieve results with a further five-fold reduc- 
tion in CPU time. 

Grid Convergence 

The effect of grid resolution in the body-normal direction 
was investigated for both the LAURA and UPS Case 2, 
two degree angle of attack solutions. The two grids used 
the same number of cells in the streamwise and circum- 
ferential directions, but had 64 and 128 cells in the body- 
normal direction, respectively. The wall cell size for the 
128 cell solutions was set to be half that of the 64 cell 
solutions, so as to maintain the same grid stretching. A 
tighter wall clustering with 64 cells was found to result 
in an over-stretched grid. 

Figures 19 and 20 plot the windside and leeside cen- 
terline heat transfer rates from the two LAURA solu- 
tions. Heat transfer rates, being a gradient of the nu- 
merical solution, are particularly sensitive to variations 
in the solution, and are thus considered a good indication 
of how well a calculation has converged. In both of these 
plots the LAURA solution can be seen to vary by 20-25 
percent in the heating between the two grids. Obviously, 
for these particular conditions LAURA is not grid con- 
verged with 64 cells in the body-normal direction. As 
seen before in Figs. 16 and 17, the 128 cell LAURA solu- 
tion agrees well with the UPS and PNSNQ3D solutions, 
so the baseline Case 2 results use the 128 cell LAURA 
solution. 

The corresponding UPS heat transfer rates for the 
two grids are presented in Fig. 21, for the windside cen- 
terline, and Fig. 22, for the leeside centerline. It is im- 
mediately apparent from these heating plots that the 
UPS solution was grid converged with the 64 cell grid. 

On both the windside and leeside there is a difference 
in heating between the grids in the immediate vicinity 
of the interface, but this is because the two solutions 
were started from the corresponding LAURA solutions, 
which displayed a significant difference in heating on 
the two grids. It is interesting to note that while the 
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Fig. 19 Case 2 body-normal grid resolution: LAURA 
windside heating. 



Fig. 20 Case 2 body-normal grid resolution: LAURA 
leeside heating. 



Fig. 21 Case 2 body-normal grid resolution: UPS 
windside heating. 



Fig. 22 Case 2 body-normal grid resolution: UPS 
leeside heating. 


starting planes for the two UPS solutions were different, 
within only five nose radii downstream the UPS solutions 
have converged, indicating that that UPS is relatively 
robust with regards to the blunt.-nose starting solution 
and grid distribution. It had been expected that LAURA 
would have been grid converged with fewer points than 
UPS, since it employs a second-order accurate numeri- 
cal scheme in all three dimensions. Figures 19-22 show 
clearly that, in fact, the opposite is true for this partic- 
ular case. 

The grid convergence of the UPS marching step size 
in the present method for Case 2 was investigated by re- 
peating the solution with a step size of 0.05 nr, five times 
larger than was used for the baseline solution. Wind- 
side and leeside centerline comparisons of heat transfer, 
Fig. 23, and surface pressure, Fig. 24, are presented for 
both step sizes. Clearly, the baseline UPS portion of 
the present method’s solution is grid converged with re- 
spect to marching step size at 0.01 nr. The cone-cylinder 
and cylinder-flare junctions are slightly better resolved 
for both the heating and surface pressure for the smaller, 
0.01 nr, step size solution, which would be expected. The 
0.05 nr step size F PS solution was obtained with 86 steps 
in 169 CPI! seconds on a Cray YMP. 

Five Degrees Angle of Attack 

A further nonequilibriunr air solution for the PNSNQ3D 
code is presented in Ref. 32 for the Case 2 configura- 
tion, but at five degrees angle of attack. Figures 25 
and 26 present the corresponding windside and leeside 
heat transfers, respectively, for the present method along 
with the results of Bhutta, who only reported a PNS so- 
lution and not a VSL solution for this configuration. 

Leeside agreement is excellent between all three so- 
lutions, with the present method matching the full-body 
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Fig. 23 UPS Case 2 marching step size convergence 
check — heating . 



Fig. 24 UPS Case 2 marching step size convergence 
check — surface pressure. 



Fig. 25 Windside heat transfers at five degrees angle 
of attack. 



Fig. 26 Leeside heat transfers at five degrees angle of 
attack. 



Fig. 27 Surface pressures for Case 2, five degrees 
angle of attack. 


LAURA solution to within three percent. The windside 
agreement is fair, though not as good as the leeside. On 
the windside centerline the UPS heating is seen to be 10 
percent higher on the cone, seven percent higher on the 
cylinder, and 18 percent higher on the flare. Also, the 
UPS and LAURA heating trends appear to be separating 
at the tail end of the body. 

The windside and leeside surface pressures from 
LAURA and UPS are shown in Fig. 27. For the sur- 
face pressure excellent agreement is seen on the windside,, 
while good agreement is seen on the leesid®;, a slightly 
different trend than for the heating. Also, the two solu- 
tions are in very good agreement on the surface pressure 
at the end of the body, as contrasted with the windside 
heating trend in Fig. 25. 

The LAURA convergence history is plotted in 
Fig. 28. Again, the early spikes in Fig. 28 result from 
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Fig. 28 LAURA convergence history for Case 2, five 
degrees angle of attack. 


grid alignments whilf the later spikes are caused by the 
multi-tasked restart procedure in LAURA. The LAURA 
solution was achieved with some difficulty for this case. 
Because of the strong expansion on the leeside at the 
cone-cylinder junction the solution had to be relaxed 
very conservatively to maintain stability. The full-body 
solution required 20 CPU hours on a Cray YMP. By con- 
trast, the LIPS solution required only 800 CPU seconds 
on the same machine, nearly two orders of magnitude 
less time. 

Concluding Remarks 

A new procedure has been implemented for the 
aerothermodynamic solution of hypersonic, chemically- 
reacting air flowhelds that combines two proven, existing 
solvers. The robustness of the thin-layer Navier-Stokes 
solver LAURA has been joined with the speed of the 
parabolized Navier-Stokes solver UPS. The class of ve- 
hicles to which the method is applicable are blunt-nosed 
configurations with slender afterbodies. The method of- 
fers the potential benefits of obtaining efficient solutions 
with second-order accuracy in the crossflow planes, while 
requiring only a fraction of the computer time and mem- 
ory that a full-body LAURA solution would require. 

Surface pressure and heat transfer results from the 
present method compare well with the baseline LAURA 
solution for the first case considered, an axisymmetric six 
degree sphere-cone at Mach 25. The downstream solu- 
tion to 100 nose radii with the present method compares 
well with the surface pressure of a viscous shock-layer 
solution, but the viscous shock-layer heating is as much 
as 20 percent lower than the present method. For the 
non-catalytic wall boundary condition it was found that 
the differences in chemistry models between LAURA and 
UPS created a stability restriction on the marching step 


size of the UPS solution, which tended to offset the de- 
crease in solution time expected with a marching scheme. 

The second case considered, a blunted multi-conic 
at Mach 28, showed good agreement between the present 
method and an all-body LAURA solution for surface 
pressures and heat transfers. Results were obtained at 
both two and five degree angles of attack. The effect 
of grid resolution was investigated in the body-normal 
direction for both UPS and LAURA, and in the stream- 
wise direction for the UPS solution. This case employed 
a fully-catalytic wall boundary condition, and did not en- 
counter any stability restriction on marching step sizes. 
A significant reduction in both computer memory and 
solution time is demonstrated with the present method 
over an all-body thin-layer Navier-Stokes solution. 

Further extensions of the present method would al- 
low for a transition from UPS back to LAURA to han- 
dle embedded separation bubbles and stagnation points. 
Also, it is anticipated that the present method might be 
used to quickly initailize more complex flowhelds, which 
would then be fully converged with a full-body LAURA 
solution. 

The present method is shown to be a fast, efficient 
procedure for obtaining aerothermodynamic predictions 
on blunted, slender vehicles at hypersonic speeds with 
reacting air howhelds. 
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